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Abstract. The angular distribution of the final electrons in the so-called long range 
mechanism of the neutrinoless double beta decay (Oi'2/3) is derived for the gen- 
eral Lorentz invariant efi'ective Lagrangian. Possible theories beyond the SM are 
classified from their effects on the angular distribution, which could be used to 
discriminate among various particle physics models inducing 0v2f3 decays. How- 
ever, additional input on the effective couplings will be required to single out the 
light Majorana-neutrino mechanism. Alternatively, measurements of the effective 
neutrino mass and angular distribution in 0v2l3 decays can be used to test the corre- 
lations among the parameters of the underlying physics models. This is illustrated 
for the left-right symmetric model, taking into account current phenomenological 
bounds. 

1 Introduction 

Neutrinoless double beta decay {Gv2(3) is forbidden in the Standard Model (SM) 
by lepton number (LN) conservation, which is a consequence of the renormal- 
izability of the SM. However, being the low energy limit of a more general 
theory, an extended version of the SM could contain nonrenormalizable terms 
(tiny to be compatible with experiments), in particular, terms that violate LN 
and allow the 0i^2/3 decay. Probable mechanisms of LN violation can include 
exchange by: Majorana neutrinos z/^/S [1-3] (one of the main candidates after 
the observation of neutrino oscillations [4]), SUSY Majorana particles [5-10], 
scalar bilinears (SBs) [11], e.g. doubly charged dileptons (the component 
of the SU{2)l triplet Higgs scalar etc.), leptoquarks (LQs) [12], right-handed 
Wu bosons [3, 13] etc. From these particles light i^s are much lighter than the 
electron and others are much heavier than the proton, giving rise to the two 
possible classes of mechanisms for the 0^2^ decay called the long range and the 
short range mechanism, respectively. For both the classes, the separation of the 
lepton physics from the nuclear physics takes place [14] which simplifies calcu- 
lations. For the first class, in contrast to the second class, the pion exchange 
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mechanism is suppressed, and in this case operators (in the effective field the- 
ory) stemming from the Hght neutrino exchange have precisely the same form 
as the leading order heavy particle exchange 01/2/3 decay operators, enabling 
a precise comparison among models [15]. According to the Schechter-Valle 
theorem [16] any mechanism of the 0i>2(3 decay produces an effective Majo- 
rana mass for the neutrino, which must therefore contribute to this decay in 
any case. These various contributions will have to be disentangled to extract 
information from the Ov2(3 decay on the characteristics of the sources of LN 
violation, in particular, on the neutrino masses and mixing. 

Despite a lack of confirmation for the claimed observation of the 0i/2f3 de- 
cay [17], the restrictions on the decay half-life [18] make it possible to get bounds 
on the parameters of the models with LN violation (see [14] for a recent dis- 
cussion). Once the Op2(3 decay has been established with good accuracy in 
the forthcoming experiments, the characteristics of this decay (half-life and the 
angular correlation between the electrons) could be combined with all the up- 
dated information from other experiments on the neutrino mixing and masses 
(neutrino oscillations, tritium beta decay [19], cosmology (WMAP [20]) etc.) 
to perform a best fit in the multidimensional space of parameters of a general 
underlying particle physics model. The fit parameters also include the masses 
and couplings of the nonstandard particles that could be involved in various 
LN- violating processes mentioned earlier. This would allow to determine the 
dominant mechanism (or a set of competing ones) of the 0^2/3 decay. 

Our aim in this paper is to examine the possibility of the determination 
of the decay mechanism from the angular correlation of the final electrons in 
this process. We restrict ourselves to the long-range mechanism and derive the 
angular distribution for the general Lorentz invariant effective Lagrangian. The 
experimental facilities that can measure the electron angular distribution in the 
Oi/2/3 decay are NEM03 [21], ELEGANT V and others [22]. We argue that 
the measurement of the angular correlation coefficient in these experiments 
would provide discrimination among the various competing scenarios of the 
Oi/2/3 decay. We illustrate this by parametrizing the angular distribution as 
dT/d cos 9 ~ 1 — KcosO (K = 1 for the hght Majorana mass scenario). Using 
the example of the left-right symmetric model [23] , we work out the correlation 
among the angular coefficient K, the mass of the right-handed Wr boson, mwR, 
and either the effective Majorana neutrino mass (m) = J2i Ueimi or the half- 
life Ti/2, taking into account the current bounds on the various parameters. It 
is shown that for values of |(m)| below 10 meV, the angular correlation between 
the electrons distinguishes the left-right symmetric models from the SM -|- light 
Majorana mass scenario. 
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2 Angular distribution for the long range mechanism of 0^2/3 decay 

2.1 General effective Lagrangian 

For the decay mediated by hght i^ms, the most general effective Lagrangian 
is the Lorentz invariant combination of the leptonic ja and the hadronic Ja 
currents of definite tensor structure and helicity [24] 

where the hadronic and leptonic currents are defined as: = uOad and 
= eOpVi; the leptonic currents contain neutrino mass eigenstates and the 
index i runs over the light eigenstates. Here and thereafter, a summation 
over the repeated indices is assumed; a,(}=V^A,S^P,TL^R iPxp = 2a'^Tp, 
a^" = I [7'', 7"^], Pp is the projector, p = L, R); the prime indicates the sum- 
mation over all Lorentz invariant contributions, except for a = (3 = V — A, and 
Uei is the PMNS mixing matrix [25]. Note that in Eq. (1) the currents have 
been scaled relative to the strength of the usual V — A interaction with Gf 
being the Fermi coupling constant. The coefficients e^^ encode new physics, 
parametrizing deviations of the Lagrangian from the standard V — A current- 
current form and mixing of the non-SM neutrinos. 

In discussing the extension of the SM for the 01/2(3 decay, Ref. [3] considered 
explicitly only nonstandard terms with 

^V-A _ g'vrr/ V + A _„T// ^V + A _x9'v^ /r,N 

^V+A,i — '^~~'^ei^ ^V-A,i — 'l^ei' ^v+A,i — '^ei ■ {■^) 

9v gv 

Implicitly, also the contributions encoded by the coefficients Cyl^ j are dis- 
cussed arising from the non-SM contribution to Uei in SU{2)l x SU{2)fi x C/(l) 
models with mirror leptons (see Ref. [3], Eq. (A. 2. 17)). Here V, U' and V are 
the 3x3 blocks of mixing matrices for non-SM neutrinos, e.g., for the usual 
SU{2)l X SU{2)b X C/(l) model V describes the lepton mixing for neutrinos 
from right-handed lepton doublets; for SU{2)l x SU{2)ii x U{1) model with 
mirror leptons [26] U' (V) describes the lepton mixing for mirror left(right)- 
handed neutrinos [3] etc. The form factors gv and g'y arc expressed through 
the mixing angles for left- and right-handed quarks. Thus, gv = cos 9c and 
g'y = e*'' cos O'q, with 6c being the Cabibbo angle, 9'(j is its right-handed mixing 
analogouc, and the CP violating phase 5 arises in these models due to both the 
mixing of right-handed quarks and the mixing of left- and right-handed gauge 
bosons (see Ref. [3], Eq. (3.1.11)). The parameters k, t], and A characterize 
the strength of nonstandard effects. Below, we give some illustrative examples 
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relating the couplings ey_^^^, ey±^j and the particle masses, couplings and 
the mixing parameters in the underlying theoretical models. 

In the R-parity- violating (RPV) SUSY accompanying the neutrino exchange' 
mechanism [5-10], SUSY particles (sleptons, squarks) are present in one of 
the two effective 4-fermion vertices. (The other vertex contains the usual Wl 
boson.) The nonzero parameters are 



^S+P 



where the index n runs over e, /i, r (1, 2, 3), and the RPV Minimal Super- 
symmetric Model (MSSM) parameters t?s depend on the couplings of the RPV 
MSSM superpotential, the masses of the squarks and the sleptons, the mixings 
among the squarks and among the sleptons. Concentrating on the dominant 
contributions ^s+Pi ^^'^ ^Tr i others are helicity-suppressed) , one can 

express Tl^q^Ln ^"^"^ ^(7)Lii '"^ follows [9] 



rjni = V sin 20f,. ( -4 ^] , (4) 

where k is the generation index, 0^^,^ and are the squark and slcpton mixing 
angles, respectively, mj^ and m^^ are the sfermion mass eigenvalues, and Ajjfe 
and X^jf. are the RPV-couplings in the superpotcntial. 

For the mechanism with LQs in one of the effective vertices [12] , the nonzero 
coefficients are 



(5) 





where 



f = Ueicii, (6) 



the parameters es(v'), Q^l^^yj, ay^yj depend on the couplings of the rcnormaliz- 
able LQ-quark-lepton interactions consistent with the SM gauge symmetry, the 
mixing parameters and the common mass scale Mg^^v) of the scalar (vector) 
LQs [27]. 
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The upper bounds on some of the parameters (6) from the Heidelberg- 
Moscow experiment were derived in Rcf. [28] using the s-wavc approximation 
for the electrons, considering nucleon recoil terms and only one nonzero pa- 
rameter e^^ in the Lagrangian at a time, see Table 1. 

Tabic 1; Upper bounds on some of parameters (6) (CL = 90%). 



y+A 


y+A 


^S+P 


'^S-P 


^" u 




6 X 10" ^ 


4 X 10"^ 


9 X 10"^ 


9 X 10"^ 


1 X 10"^ 


6 X 10"'^ 



The coefficients e^^ entering the Lagrangian (1) can be expressed as 

J - 



(7) 



where Uei are mixing parameters for non-SM neutrinos (see, e.g., Eq. (2)). As 
this Lagrangian describes also ordinary /3-decays (without LN violation), the 
coefficients are constrained by the existing data on precision measurements 
in allowed nuclear beta decays, including neutron decay [29]. For example, 
from these data we obtain the conservative bound 



< 7x10-2. 



(8) 



Prom Eqs. (6), (7), (8) and Table 1 we can assume that the nonstandard mixing 
is small: 

'UeiUei <10-'. (9) 



2.2 Approximations and electron angular distribution 

We have calculated only the leading order in the Fermi constant and the leading 
contribution of the parameters to the decay width using the approximation 
of relativistic electrons and non-relativistic nuclcons. Following Rcf. [2] we de- 
scribe the outgoing electrons by plane waves approximately taking into account 
the effect of the nuclear Coulomb field by the Fermi factors F{es) [2,3] with the 
s-th electron energy Eg- The non-relativistic structure of the nucleon currents 
in the impulse approximation is taken from Ref. [30]. We neglect the nucleon 
recoil terms, because the accurate calculation of the corrections should also in- 
clude simultaneously a precise calculation of the effect of the nuclear Coulomb 
field, requiring the technique of the spherical waves for the electrons. Note that 
in Ref. [30] the recoil terms due to the weak magnetism were calculated in the 
model with only V^A currents (the recoil terms due to the pseudoscalar form 
factor were not taken into account) and it was shown that in the interaction 
proportional to e^l^ the recoil effect dominates over other contributions in 
the 0+ 0+ transition. However, our numerical calculations carried out in 
section 3 are restricted to the coefficients Uei (i-e., the SM -|- light Majorana 
us scenario) and Cyt^ j (light Majorana us in the left-right symmetric models) 
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entering the Lagrangian (1), and in both cases the nucleon recoil effects are not 
dominant. 

We obtain the differential width in cos 9, where 9 is the angle between the 
electron momenta in the rest frame of the parent nucleus in the 0'^{A, Z) — > 
0"'"(A, Z + 2)e~e~ transitions, 

C\MGT\^I[{a + h){l-k COS 6)] , 



dcos9 

I[x] = / dsielelFiei)F{e2)x , (10) 
where £2 = A — £1 and A is the energy release in the process. The constant 

contains the electron mass rUe and the nuclear radius Rq, included in the def- 
inition of C so that the a and b functions and the neutrino potentials are 
dimensionless. The Gamow-Teller nuclear matrix element, 

Mgt = (0+||^/i(r„6,a;)tT„ • tT6r»T|||0+), (12) 

contains the neutrino potential h{r,u}) = Ro<j)o/r with (/>o = e*"*", r = Tab is 
the distance between the nucleons a and h, and w is the average energy of 
the neutrino. The operator t° = (ti + zr2)°/2 converts the a-th neutron into 
the a-th proton; |0^) ((0^ |) is the initial (final) nuclear state. The angular 
correlation coefficient in Eq. (10) is 

fc = ^— -1<A;<1. (13) 

The expressions for a and b for different choices of e^, considered only one 

at a time, are shown in Table 2. In this table £12 = £1 — £2 and (m) is 
the effective Majorana mass. The form factors gv{(f), dAiq^), Fs^\q^)j and 
T['^\q^) describe the following nucleon matrix elements [31] 

{P{k')\ud\N{k)) = Ff\q^){[^{k')T+i,{k), (14) 
{P{k')\u2^^'PLMN(k)) = i>(k')^^ [gv{.q^) T 9A{q^)l^] r+^P{k), (15) 

{P{k')\u2a^''PL,Rd\N{k)) = i,{k') 



T+i}{k), 
(16) 

where 
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is a nucleon isodoublet. We neglect the dipole dependence of the form factors 

on the momentum transfer q = k' — k and omit the zero argument of the form 
factors. Notations similar to the ones in Ref. [3] are used in Table 2. 

Tabic 2: Expressions for a and h in Eqs. (10) and (13) for tli(> stattxl choic(> of e£. 



( 


(1 


b 


V-A 
^V-A 


\[{m)/me + 2 Ue.el.Zi ,irn,, / nwMl - Xf)\'' 





y-A 

^V+A 


1(1 - XF){m)/me - 2 C/e^4;l,K/me)(l + Xf)\' 





^S+P 


ao + i(A/me)2(Ff /5y)2|e|i^x'Fp 







ao + i (Tf V5A)^(A/me)2|e?«(x'GT + ^x'tW 







ao 





e 


b 


a 


y^A 

^V-A 


||e;;i^R(ei2/me)2|x2+P + UA/mer\x'p\'] 


ao 


y+A 


U^i2/merKX^X2-\'' 


«o 




^(^'"/«/r)-| Uc,(stpA"''/>'>^ )Xj-1- 


flo 


4r 


52{Ti'^/gAr\j:iUeielU^i/me)? 


ao 



Thus, 



X2± = XGTo; iXFw - gXlT' Xl± = x'gT - ^X't ^ ^X'f' (18) 



Xf= \ — 

.9 A 



Mf 



XP = 1 



M 



Mgt \9aJ Mgt Mgt 

with Fermi Mp, pseudoscalar Mp, and tensor Mt nuclear matrix elements: 

Mf = (0+||^Mr.fe,u;)r>|||0+), (20) 



Mp = {0+\\J2Krat„co) {{aa - ab) ■ [n x n+]}T»T^||0,+) 

1 



(21) 



MT = {Oj\\^hirab,i^) 



a a ■ no-fo • n - -cTa • 



T;rtl|0+), (22) 



with 



r = rn, R = i?n+, (23) 

where r = r,,;, (R = R„6) is the difference (half sum) of radius-vectors of 
the nucleons a and 6; n and n_|_ are unit vectors. The prime and the index 
(J imply that the matrix element in the numerator instead of h contains the 
neutrino potential h! = h + LoR^hi or h^^ — h — ujR{)hi, respectively, with 
hi = —d(j)o/d{u)r). The quantity a for all zero is called ao in Table 1, and is 
defined as: 

ao = |(l-XF)(m)/me|'. (24) 



For the V ^ A part of the Lagrangian (1), our result agrees with Ref. [3] for 

the relativistic electrons (rrie/A 0) that weakly interact with the nucleus 
{aZ ^0), if the recoil and P-wave effects are not taken into account. 

If the effects of all the interactions beyond the SM extended by the ums, 
which we call the "nonstandard" effects, arc zero (i.e., all =0), then k = 1 
and the distribution (10) is proportional to 1 — cos9. The angular coefficient 
deviates from 1 only for the cases b ^0 irrespective of the value of a. Therefore, 
the presence of the "nonstandard" first set of parameters in Table 2, 



^S+P Tn 



.V-A 



and e^" does not change the form of the angular distribution, 
but the presence of the second set (see lower part of Table 2), e^^^, e|~p, e^^, 
e^^ and e^^ does change this distribution. Thus, experimentally establishing 
fc 7^ 1 would signal the presence of beyond-the-SM contribution in the Oj/2/3 
decay. The converse is not true; namely establishing k = 1 experimentally 
will not single out the SM + um^ as the only mechanism of the decay 
and one would require additional input/constraints on the parameters of the 
underlying theory with their coefficients listed in the upper part of Table 2. 
The coefficient k and the set {e} of nonzero e^s that change the 1 — cos^ form 
of the distribution for the SM plus i^ms are given in Table 3 (the lower two 
entries). They correspond to the following extensions of the SM: i/ms plus RPV 
SUSY [9], VMS plus right-handed currents (RC) (connected with right-handed 
W bosons [3] or LQs [12]). 

Table 3: The angular correlation coofRcicnt k for various SM extensions. 



SM extension 




k 






1 


i/M+RPV SUSY 


S-P 
'^S+P 


1 -fc< 1 




y+A 


-f < fc < 1 



Among the models that we have listed in Table 3, the coefficient e|_)_p for 
the vm-\- RPV SUSY case is helicity suppressed (and eflp = 0), and hence 
the angular coefficient A: ~ 1 for this model. (It is a mathematical challenge to 
come up with a model which lifts this chiral suppression). Among the realistic 
models we have discussed, only the model called vm-\- RC can essentially change 
the angular coefficient k from being 1. Left-right symmetric models belong to 
this class and we have studied these models in detail in section 3, where the 
correlation among the parameters K (see Eq. (34) below), mwa and |(m)| (or 
T1/2) is worked out. 

For a quantitative understanding of the influence of the decay mechanism 
on the angular distribution, it is necessary to take into account the effect of 
the nuclear Coulomb field on the electrons in the terms with efj n and e!!^^ " 

O^P J-L.R 

and the nucleon recoil terms. Note that the calculation of the above-mentioned 
corrections for the long range mechanism of the 02/2/3 decay will change inessen- 
tially the values of a and h for Cyl^, Cy+^A or Cy^^, varied one at a time, and 
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hence the results of the next section. 



3 Electron eingular correlation in left-right symmetric models 

The experimental bounds on the axe connected with the masses of new parti- 
cles, their mixing angles, and other parameters specific to particular extensions 
of the SM [2, 3, 7, 9, 11, 12]. To illustrate the kind of correlations that the mea- 
surements of \{m) \ and the angular correlation coefficient k in the 0^2(3 decay 
would imply, we work out the case of the left-right symmetric models [23]. In 
the model SU{2)l x S'[/(2)ij x U{1) the parameter A (see Eq. (2)) is expressed 
through the masses and of the left- and right-handed W bosons [3]: 

X={mwL/mwRf, (25) 

under the condition 

< mwR- (26) 

Eqs. (2) and (6) yield the relation 

e^li = X^UeiVei. (27) 

9v 

To reduce the number of free parameters, we assume the equality of the form 
factors of the left- and right-handed hadronic currents: 

9v = g'v (28) 

The small masses of the observable vs are likely described by the seesaw formula 
that in the simplest case gives 

m^r^ml,/MB., MR:^mD, (29) 

with the Dirac mass scale (for the charged leptons and the light quarks 
mo }^ 1 MeV) and the mass scale Mr of right vms (in the majority of theories 
Mr > 1 TeV). In the left-right symmetric models these scales arise usually 
from the two scales of the vacuum expectation values of Higgs multiplets [23]. 
In the seesaw mechanism, the values of the mixing parameters Ki (for i num- 
bering light mass states) have the same order of magnitude as mo /Mr. In our 
discussion we use two rather conservative values (compare with Eq. (9)) 

e = 10-^ 5 X 10"^ (30) 

for the mixing parameter 

e = \UeiVei\. (31) 
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We recall that here the summation index i runs only over the light neutrino 

mass eigenstates (the summation over the total mass spectrum including also 
heavy states gives strictly zero due to the orthogonality condition [3]). 
Prom Eqs. (25), (27), (28) and (31) we have 

mwR = mw^ (e/ |ey^l|)^^^ • (32) 

Using Eq. (26) we note the approximate equality of itiwl and the mass of the 
observed charged gauge boson Wi (mvyi=80.4 GeV [4]). 

We now turn to work out the relations among the angular correlation coeffi- 
cient fc, the right-handed W^-boson mass and the neutrino effective mass 
|(m)|. To this end, we note that the differential width for the only nonzero 
nonstandard parameter e^^^ can be obtained by comparing Eqs. (10), (13), 
(24) and the corresponding expression for b from Table 2 with the more precise 
result of Ref. [3] , yielding: 

= \Mgt\'^-^{A + B){1 -Kcose), (33) 

which is correct for the spherical waves of electrons distorted by the Coulomb 
field of the nucleus in the limit of small rUe/A. Here, 

^=J^' A = aoGoi, B = \e^Xix2-\^ Go2, (34) 

with the usual phase space factors Goi {i — 1,2) defined in Ref. [3]. Note that 
the angular coefficient K entering in Eq. (33) differs from the coefficient k, 
entering Eqs. (10) and (13), as is a function of the electron energy £i and K 
is obtained by integrating over the energy spectrum in Eq. (10). The ag and X2- 
in Eq. (34) contain instead of (po the neutrino potential (f> from Ref. [3] . In the 
numerical calculation we have used xf = 0.274, X2- — 0.551, Mqt = 1.846, 
obtained in a QRPA model with p-n pairing for ''^Ge, with Gqi = 7.928 x 10^^^ 
and Go2 = 12.96 x 10~^^ {y^~^) [32]. For a recent discussion on uncertainties 
in 01^2/3 decay nuclear matrix elements see Ref. [33]. 
Using Eqs. (34), (24), and (32), we have 

y + 1' Go2\ X2-e rUe J \mwL J 

The correlation among K, niwa and |(m)| is shown in Fig. 1 for e = 10~^ 
and in Fig. 2 for c = 5 x 10^'''. We consider the values of \{m)\, starting from 
the current upper bound from the Heidelberg-Moscow experiment, taken as 
I (to) I < 0.3 eV, up to I (to) I = 0.001 eV, covering most scenarios of neutrino 
mass hierarchies and mixing angles (see Ref. [34] for a recent discussion and 
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update). Concerning the existing bounds on mwR, we note that from Eqs. (2), 

(7) and (28) one obtains e^^^ = A. With this, Eq. (25) and the constraint 

(8) derived from [29] yield niwR > 300 GeV. This bound is weaker than the 
one mwfi > 715 GeV, obtained from the electroweak fits [4]. There is still 
a more stringent bound mwn > 1-2 TeV, obtained in Ref. [35] for the decay 
mediated by heavy Majorana neutrinos using arguments based on the vacuum 
stability [5], but it requires additional theory input. We assume tuwr > 1 TeV 
in all our figures. 

Using (32), (33), (34) and the relation T1/2 = ln2/r we get the correlation 
among tuwr and the measurable 01/2/3 decay parameters, namely the half-life 
Ti/2 and the angular coefficient K: 

K=l- 2Go2{\MGT\X2-ef{mwJmwnfT^/2- (36) 

The correlation among K, uiwr and T1/2 is shown in Fig. 3 for e = 10~^ and 
in Fig. 4 for e = 5 X 10^'''. 

Figs. 1-4 are the principal numerical results of this paper. They show 
that depending on the values of |(m)| (or T1/2) and ttiwr, all values of the 
angular coefficient K are allowed. For example, Fig. 1 shows that it is possible 
to describe the angular distribution close to 1 + cos 6 {K > —1) by the (long 
range) mechanism with the right-handed boson Wr with the mass about 1 TeV 
for I (to) I ^ 1 mcV. 

For illustration, in Fig. 5 we plot the differential width (33) vs. cos9 for 
a set of values of |(to)| and itiwr, assuming e = 10~^. It is seen that the 
sensitivity of the electron angular distribution to the right-handed VK-boson 
mass mwR increases with decreasing values of the effective Majorana neutrino 
mass I (to) I, as can be seen from Fig. 5 (right), where this distribution is shown 
for I (to) I = 1 meV, 3 meV and 5 meV. 
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Figure 1: Left: Dependence of the angular correlation coefficient K on the right-handed W- 
boson mass and the value of the neutrino effective mass for the 0;/2/3 decay of 

'^^Ge. The mixing parameter e = 10~®. Right: The same as the left figure but for smaller 

values of |{m)|. 
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Figure 3: Left: Dependence of the angular correlation coefficient K on the right-handed W- 
boson mass and the half-life Tiy2 for the 0j/2j8 decay of ''^Ge. The mixing parameter 

e = 10~^. Right: The same as the left figure but for larger values of 
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Figure 5: Left: Dependence of the differential width for the 0u2f3 decay of '°Ge on cos0 for 
a fixed value of e = 10~^ and \{m)\ = 10,15,20,30 meV. The dashed, dotted and straight 
lines correspond to m-^^ = 1, 1.2,oo TeV, respectively (the latter is the conventional case 
of the light Majorana neutrino exchange mechanism). Right: The same as the left figure but 
for smaller values of = 1,3,5 meV. In addition, the dash-dotted lines correspond to 

mw„ = 1.5 TeV. 
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